Magnetic Susceptibility of Collinear and Noncollinear Heisenberg Antiferromagnets 
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Predictions of the anisotropic magnetic susceptibility % below the antiferromagnetic (AFM) order- 
ing temperatures Tn of local moment Heisenberg AFMs have been made previously using molecular 
field theory (MFT) but are very limited in their applicability. Here a MFT calculation of %(T < Tn) 
is presented for a wide variety of collinear and noncollinear Heisenberg AFMs containing identical 
crystallographically equivalent spins without recourse to magnetic sublattices. The results are ex- 
pressed in terms of directly measurable experimental parameters and are fitted with no adjustable 
parameters to experimental x(T < Tn) data from the literature for several collinear and noncollinear 
AFMs. The influence of spin correlations and fluctuations beyond MFT is quantified by the devi- 
ation of the theory from the data. The origin of the universal x(T < Tn) observed for triangular 
lattice AFMs exhibiting coplanar noncollinear 120° AFM ordering is clarified. 



Introduction. Magnetic susceptibility x measure- 
ments versus temperature T have been used for a cen- 
tury to obtain important information about the magnetic 
properties of materials. The Weiss molecular field theory 
(MFT) has been instrumental in interpreting the x(^) 
data in the paramagnetic state above the long-range mag- 
netic ordering temperature Tn of local magnetic moment 
antiferromagnets 1,2 (AFMs) via the Curie- Weiss (CW) 
law x — T-e •> m which the magnitude of the local mo- 
ments is contained in the Curie constant C and the na- 
ture and strengths of their interactions in the Weiss tem- 
perature P . MFT has also been used extensively for 
comparisons with experimental data of its predictions for 
the ordered magnetic moment and magnetic heat capac- 
ity versus T in the ordered state of AFMs at T < T N . 
Thus MFT is a primary tool to identify important char- 
acteristics of local moment AFMs. 

In contrast, very few comparisons have been made of 
experimental anisotropic x(T < Tn) data for AFMs with 
the predictions of MFT even for collinear AFMs where 
the ordered moments jli are aligned along the same easy 
axis. 1 3 Here we provide simple MFT expressions to fit 
experimental x(T < Tn) data for ordered AFMs con- 
taining identical crystallographically equivalent spins in- 
teracting by Heisenberg exchange for arbitrary sets of 
exchange constants. The theory treats collinear and pla- 
nar noncollinear AFM structures on the same footing 
without the use of magnetic sublattices. The results are 
expressed in terms of independent experimentally mea- 
surable quantities and are used to fit with no adjustable 
parameters representative experimental x(T < Tn) data 
from the literature for several collinear and noncollinear 
AFMs. The fits can quantify the influence of spin cor- 
relations and fluctuations beyond MFT on x(T < Tn), 
and can also help to elucidate the AFM structures and 
exchange interactions if these are uncertain or unknown. 

Using MFT, Van Vleck calculated in 1941 the 
anisotropic x(T < Tn) for magnetic fields H applied 
parallel (x\\) and perpendicular (x_l) to the easy axis of 
collinear AFMs with only nearest-neighbor Heisenberg 
interactions between spins on two distinct interpenetrat- 
ing "bipartite" sublattices. 4 Yoshimori carried out MFT 



Proper Screw Helix 
Side View Top View 




X 



FIG. 1: (Color online) Yoshimori's "proper screw helix" non- 
collinear AFM structure that he proposed for MnCb- 5 Each 
magnetic moment represents a planar array of ferromagnet- 
ically aligned magnetic moments that lie in the xy-pl&ne, 
where the planes are separated along the helix axis (z-axis) 
by distance d. The plane of the magnetic moments is perpen- 
dicular to the helix wave vector k along the z-axis. The net 
interplanar exchange interactions J z ± and J Z 2 in the generic 
Jo-Jzi-Jz2 model are indicated. A "cycloidal helix" AFM 
structure 5 occurs when k is in the xy-p\&ne of the magnetic 
moments. 



calculations of x(T < T^) in 1959 for the special case 
of a planar noncollinear AFM "proper screw helix" mag- 
netic structure that he proposed for Mn02, 5 ' 6 as shown 
schematically in Fig. 1. These MFTs are very restricted 
in their applicability and have been rarely used to fit ex- 
perimental x(T < Tn) data over the past five decades. 

Theory. Here we consider identical crystallograph- 
ically equivalent spins interacting by Heisenberg ex- 
change, with no anisotropy present except that due 
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to an infinitesimal H. The part Ei of the aver- 
age energy of the system that is associated with in- 
teractions of jli with H and with its neighbors jij 
is Ei = Wj^fc ' Ji ofij - Mi ■ H, where is the 
Heisenberg exchange coupling between ordered magnetic 
moments jli and jij. Using MFT, 1 one obtains the 

CW law for T > T N , where C = , # P = 

-^^EjJij and T N = -?<m ZjJijMS^ N is 
the number of spins, g is the ^-factor, /ip is the Bohr 
magneton, S is the spin, &b is Boltzmann's constant and 
<j)ji are the angles between jli and its neighbors jij in the 
AFM-ordered state. We rewrite the CW law for T > T N 
in dimensionless form as 
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J2j Jij cos 4>ji' 
(i) 

Below Tn, the x with H perpendicular to the ordered 
moment axis or plane for collinear or planar noncollinar 
AFMs, respectively, is given in general by MFT as 1 



X±(T<T N )T N x(T N )T N 1 
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(f<l) (2) 



For collinear AFMs, a field applied below Tn along 
the easy axis just changes the magnitude of an ordered 
moment without rotating it and in MFT we obtain 



X\\(t)TN 
C 

where B' s (y ) = 
function, 1 yo = 



T*(t) = 



(S + l)t 



(3) 



dB s (y)/dy\ y=yo , B s (y) is the Brillouin 
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magnitude of the ordered moment in zero field p,o(t) is 
calculated numerically from /2o = -Bs(yo)- 1 From Eqs. (2) 
and (3) one obtains 



Msat = 9 Sub, and the 
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(4) 



By Taylor expanding B' s (yo) = (S + l)/3 for 7/0 ~^ 0? 
one obtains r*(£ -4- 1) = 1 and X|I ^ 7 ^"^ N ^ ) = 1, as re- 
quired. For T ->• 0, B' s (yo) ->• 0, r* ->> 00 and X|| ^ 0. 
The parameters in Eq. (4) required to fit experimental 
Xy(T < Tn) data are just /, 5 and x(Tn), which can 
usually be easily independently determined from exper- 
iment or estimated. Setting / = — 1 in Eq. (3) repro- 
duces Van Vleck's 1941 prediction for the special case 
of bipartite collinear AFMs with only nearest-neighbor 
interactions. 4 

For planar noncollinear AFMs, one must take into ac- 
count via MFT the field-induced changes in both the 
magnitudes and directions of the ordered moments to 
first order in and we then obtain the in-plane (xy) 
susceptibility 



Xxy(T < T N ) (i +r * + 2 / + 4g*)(l-/) 

X(T N ) " 2 [(r* + £*)(! + B*) - (f + £*) 2 ] : 
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where 



Ej J ii cos2 0ji 
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Using T*(t -> 1) = 1, Eq. (5) gives ^gff^ = 1, 
irrespective of the value of £?*, as required, whereas 
lim^o B' s (yo) — > and r* — > oo yield from Eq. (5) 



Xxy(T = 0) 

x(?n) 
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2(1 + 5*)' 



(7) 



The parameter 5* is the only new parameter specifi- 
cally associated with noncollinear AFMs, is not generally 
directly measurable, but can be evaluated if the AFM 
structure and an exchange interaction model are avail- 
able. Alternatively, it can be used as a fitting parameter 
to provide such information. 

On the other hand, the value of B* can be experi- 
mentally determined within a minimal generic Jo-J z i-Jz2 
model 6 for helical/cycloidal AFM structures as in Fig. 1 
on any Bravais spin lattice. In this model, one sums the 
exchange interactions of a given magnetic moment with 
all other moments in the same ferromagnetically-aligned 
layer perpendicular to the helical/cycloidal wave vector 
k and calls that sum Jo, and similarly for nearest- and 
next-nearest-layer interactions J z \ and J^2, respectively, 
as indicated in Fig. 1. The same theory is applicable to 
isolated linear chains where Jo = 0. Then the k of the he- 
lix/cycloid is obtained in terms of the exchange constants 
by minimizing the exchange energy to be 5,6 



cos(kd) 
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where k = |k| and d is the distance between layers, kd 
is the turn angle between adjacent moments along the 
helix/cycloid axis (Fig. 1) and is experimentally measur- 
able by magnetic x-ray or neutron diffraction techniques. 
Using Eq. (8) one can express B* in Eq. (6) as 



B* 



2(1 - /) cos(fcd)[l + cos(fcd)] - /. 



(9) 



Using Eq. (9), one can now write Xxy(T < Tn)/x(^n) 
in Eq. (5) completely in terms of independently measur- 
able quantities. Furthermore, using Eqs. (7) and (9) one 
obtains 



Xxy(T = 0) 

xCTn) 



1 



2 [l + 2 cos(ftd) + 2 cos 2 (kd)] 



(10) 



The expression for Xxy (T = 0)/x(T N ) obtained in 1959 
by Yoshimori 5 for the special case of the helix in Mn02 
is consistent with the general result (10). 

Using Eq. (10), the Xxy(T = 0)/x(T N ) is plotted ver- 
sus kd in Fig. 2(a). The predicted behavior has a surpris- 
ing nonmonotonic dependence on kd with a maximum at 
kd = 27r/3 with a value of unity. Using Eqs. (5) and (9), 
Xxy(T < Tn) and its dependences on kd and / are shown 
in Fig. 2(b), where Xxy is seen to be strongly dependent 
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FIG. 2: (Color online) (a) Xxy(T = 0) versus kd [Eq. (10)] 
and (b) \xy versus T and / [Eqs. (5), (9)] for helical/cycloidal 
AFMs with S = 7/2 within the generic J -J z i-Jz2 MFT 
model. 



on T and / except for kd = 2tt/3 = 120° for which it is 
independent of T and /. One can prove that this result 
for kd = 27r/3 is obtained within the Jq-J z \-Jz2 model 
for any value of S. Then using Eq. (2), our MFT makes 
the remarkable universal prediction for helical/cycloidal 
120° AFM ordering that x(T < T N )/x(T N ) is isotropic 
and independent of 5, / and T for T < Tn. The same 
result is obtained for other AFMs with 120° ordering 
and therefore a helical/cycloidal AFM structure is not 
required [see also Fig. 5(a) below]. 

If only the six nearest-neighbor interactions J occur 
in a single triangular lattice layer exhibiting 120° or- 
dering in MFT, one obtains from Eqs. (1) and (2) that 
X(T = 0)/(Ng 2 fil) = 1/(9 J), independent of S. For the 
classical (S —> oo) isolated triangular layer Heisenberg 
AFM, one obtains the same value. 7,8 Classical Monte 
Carlo simulations for a triangular spin lattice layer in- 
dicate that x is isotropic and also nearly independent 
of T at low T. 9 Our MFT result for kd = 2tt/3 thus 
significantly extends the previous calculations for single 
classical triangular lattice layers to finite quantum spins 
S and long-range AFM ordering of coupled layers. 

Fits of Experimental Data. As shown in Eq. (2), x± 
is independent of T below Tn with the value %(Tn), so 
no explicit fitting of experimental data is required. 

We first present fits by Eq. (4) of Xy(T) data for 
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FIG. 3: (Color online) Anisotropic x(T) of single crystals of 

(a) GdNiGes (Ref. 10) and (b) MnF 2 . 12,13 In (a), the easy-axis 
is the a-axis (open circles), whereas in (b) it is the c-axis. In 

(b) , the P value was taken from Ref. 14. The corresponding 
X\\ (T) data are fitted by the MFT prediction in Eq. (4) (solid 
red curves) with no adjustable parameters. 



the collinear AFMs GdNiGe3, an orthorhombic com- 
pound containing nonmagnetic Ni atoms and Gd +3 spins 
S = 7/2, 10 and MnF2 with the primitive tetragonal ru- 
tile structure containing Mn +2 spins S = 5/2. 11 The 
anisotropic x(T) data at low T for single crystals of 
GdNiGe 3 (Ref. 10) and MnF 2 (Refs. 12,13) and the cor- 
responding fits of the X\\(T < Tn) data by Eq. (4) with 
no adjustable parameters are shown in Fig. 3. The fit to 
the X|i(T < Tn) a-axis data of GdNiGe3 with S = 7/2 
is better than the fit to the corresponding c-axis data of 
MnF 2 with S = 5/2. This comparison agrees with ex- 
pectation, because MFT does not include the influence of 
quantum spin fluctuations which increase as S decreases. 
This suggests that a comparison of such MFT fits with 
experimental data is a quantitative diagnostic for the oc- 
currence at T < Tn of spin fluctuations and correlations 
beyond MFT. 

As an example of a noncollinear planar AFM, primitive 
tetragonal GdE>4 consists of crystallographically equiv- 
alent Gd spins 7/2 with the AFM structure shown in 
Fig. 4(a) and with the ordered moments oriented in 
the [110] and equivalent directions. 15 The magnetic and 
chemical unit cells are the same. Anisotropic x(T) 
data at low T are shown in Fig. 4(b). 16 The fit of the 
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FIG. 4: (Color online) (a) AFM structure of GdB 4 . 17 (b) 
Magnetic susceptibility % versus temperature T for a single 
crystal of GdB 4 , 16 together with the fit of the Xab(T < T N ) 
data by Eq. (5) using S = 7/2 and the experimentally deter- 
mined parameters / and B* in the figure. Thus there are no 
adjustable parameters in the fit. 
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FIG. 5: (Color online) Anisotropic susceptibilities Xab(T) for 
the triangular lattice AFMs (a) hexagonal YMn0 3 (Ref. 18) 
with Mn +3 spins S — 2 and kd = 120° ordering in the ab- 
plane 17 and (b) RbCuCl 3 with Cu +2 spins 5=1/2 (Ref. 23) 
that has a slightly distorted triangular spin lattice with kd — 
108° cycloidal ordering in the hexagonal afr-plane. 24 The x(T) 
data in (b) are corrected for the T-independent diamagnetic 
core and paramagnetic Van Vleck contributions. 23 



Xab(T < Tn) data by Eq. (5) with no adjustable pa- 
rameters is shown by the solid blue curve using param- 
eters in the figure. The value of B* was estimated 
from Eq. (7) and the experimental values 16 of / and 
Xab(T — > 0)/x(Tn). The relationship between the fit and 
data is similar to that for GdNiGe3 in Fig. 3(a). 

We now test our universal prediction for noncollinear 
120° AFM structures that x(T < T N ) is isotropic and 
independent of /, S and T for < T < Tn with the 
value x(Tn), which does not require explicit fits. The 
hexagonal compound a-YMn03 contains a triangular lat- 
tice of crystallographically equivalent Mn +3 spins S = 2 
and exhibits 120° coplanar ordering in the afr-plane. 17 
As in GdE>4, the magnetic and chemical unit cells are 
the same. Anisotropic x(T) data for this compound are 
shown in Fig. 5(a). 18 The x( T < t n) data parallel and 
perpendicular to the afr-plane are nearly isotropic and 
independent of T. Similar x(T < Tn) results have been 
obtained for many triangular lattice AFMs with 120° he- 
lical or cycloidal ordering, such as the 5 = 3/2 com- 
pounds LiCr0 2 , 19 VF 2 and VBr 2 . 20 22 Our MFT predic- 
tion is even strongly confirmed by the x(T < Tn) data 23 
in Fig. 5(b) for the slightly monoclinically distorted tri- 



angular spin lattice in RbCuCl3 containing highly quan- 
tum Cu +2 spins- 1/2 exhibiting cycloidal AFM order- 
ing within the hexagonal afr-plane. 24 The cycloid axis 
is in the hexagonal [110] direction with a turn angle 
kd = 108°, 24 close to the undistorted triangular lattice 
value of 120°. The reason that the MFT prediction is 
accurate even for S = 1/2 deserves further investigation. 

In summary, a generic molecular field theory of the 
anisotropic x(T < Tn) was formulated for local moment 
Heisenberg AFMs that is widely applicable to collinear 
and planar noncollinear AFM structures. The com- 
parisons of our results with experimental anisotropic 
x(T < Tn) data for single crystals in Figs. 3-5 with 
no adjustable parameters demonstrate that such anal- 
yses constitute a powerful probe of the AFM structure 
and spin interactions. Our results will also be useful for 
analyzing x(T < ^n) data for polycrystalline samples. 
An important avenue for future research is to further 
study the applicability, accuracy and limitations of our 
MFT predictions. The present work is a stepping stone 
for additional MFT calculations of x(T < Tn) that could 
include various types of anisotropics. 
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